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Abstract 

We study the vertical extent of propeller structures in Saturn's rings. Our focus 
lies on the gap region of the propeller and on non-inclined propeller moonlets. 
In order to describe the vertical structure of propellers we extend the model of 



Spahn and Sremcevic ( 2000 ) to include the vertical direction. We find that the 



gravitational interaction of ring particles with the non-inclined moonlet does 
not induce considerable vertical excursions of ring particles, but causes a con- 
siderable thermal motion in the ring plane. We expect ring particle collisions 
to partly convert the lateral induced thermal motion into vertical excursions 
of ring particles. For the gap region of the propeller, we calculate gap aver- 
aged propeller heights on the order of 0.7 Hill radii, which is of the order of 
the moonlet radius. In our model the propeller height decreases exponentially 
until viscous heating and collisional cooling balance. We estimate Hill radii of 
370m and 615m for the propellers Earhart and Bleriot. Our model predicts 
about 120km for the azimuthal extent of the Earhart propeller at Saturn's 2009 
equinox, being consistent with values determined from Cassini images. 
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1. Introduction 



Planetary rings are one of the most remarkable and beautiful cosmic struc- 



tures. They are considered to be natural dynamical laboratories (Burns and 



Cuzzi 2006), exemplifying the physics of cosmic disks, such as accretion disks 



or galactic disks, which are much larger and much farther away from Earth. An 
exciting example is the presence of small moons embedded in Saturn's rings, 
henceforth called moonlets, which have their analog in planetary embryos orbit- 



ing within a protoplanetary disk (Artymowicz 2006 Papaloizou et al. 2007). 

The fact that the resolution of even Cassini 's cameras is too low to image 
these moonlets directly, brought up the idea of investigating moonlet induced 



putative structures in the rings (Spahn 1987 Spahn and Wiebicke 1989), with 
the hope that these features could be captured by the spacecraft cameras or 



instruments. This then led to predictions of the propeller structures (Spahn 



and Sremcevic 2000 Sremcevic, Spahn, and Duschl 2002) which are carved 



in the rings by the moonlet. Subsequent numerical particle experiments (Seifi 



et al. 2005 Sremcevic et al. 2007 Lewis and Stewart 20091 completed the 



fingerprint of such gravitational perturbers and confirmed the spatial scaling of 
the propeller structure. Depending on its size, an embedded ring-moon either 
creates a propeller (sizes below 1 km) or, alternatively, opens up a circumfer- 
ential gap (for sizes > 1 km, e.g. the ring-moons Pan and Daphnis). Both 
structures, propeller and circumferential gap, are decorated with density wakes, 
completing the structural picture. Up to this stage, all these density features 
have been assumed to occur only in the ring plane, a vertical stratification of 
moonlet induced structures has not seemed to be of importance. 



More than 150 propeller moonlets have now been detected (Tiscareno et al. 



2006| 2008) and among them a few which are large enough to allow Cassini 's 



cameras to take several snapshots of their propellers at different times, con- 
firming in this way their orbital motion. Those moonlets were nicknamed after 



famous aviators, e.g.: Bleriot, Kingsford Smith, Earhart (Tiscareno et al. 2010). 
In the summer of 2009, at Saturn's equinox (the Sunset at Saturn's rings), 
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the perfect opportunity arose to detect any vertical structure deviating from 
the mean ring plane by observing shadows cast on the rings. At this time the 
density structures around the largest propeller moonlets, as well as those around 
the ring-moon Daphnis, created prominent shadows. These can be assigned to 
the wakes and in the case of the propeller moonlets also to excited regions of 
the propeller, where the moonlet induces a partial gap. The shadows were much 
longer than the moon's size itself. The lateral extent of the shadows allows to 
conclude that moonlet induced vertical excursions of ring particles can be in the 
range of several kilometers in the case of Daphnis or several hundred meters in 
the case of the large propeller moonlets. These very facts directly indicate the 
necessity to investigate the vertical stratification of moonlet induced structures, 
which has not been the focus of former models of the moonlet 's fingerprint. 
In this work we will study the vertical extent of moonlet induced propeller 
stuctures, focusing on the gap region of the propeller. 

The paper is organized as follows: In Section[2]the extended propeller model 



is presented. In Section 2.2 the mass flow through the scattering region is cal- 
culated by a probabilistic approach, and values of the moonlet induced thermal 
velocities arc determined, which are later used as initial conditions for the hydr- 



dynamical equations. Section 2.3 gives the hydrodynamical balance equations, 
which we use to model the diffusion of mass into the induced gap and the relax- 
ation of the ring temperature. In Section [3] we calculate the height of the ring 
in the gap region of a propeller. Section [4] discusses made assumptions and the 
application of our results to propeller features in Saturn's rings. Our results are 
summarized in Section [H 

2. Extended model of gravitational scattering 

2.1. The scattering region 

The first step in the formulation of our model, is to divide the planetary 
ring, composed of granular grains and one moonlet, into two regions: 

(i) the scattering region, 
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(ii) the rest of the ring. 

In this work we consider moonlets on circular, i.e. on non-inclined and 
non-eccentric, orbits. The scattering region is the small area (volume) around 
the moonlet where the majority of the trajectory changes, due to the moon- 
let's gravity, take place. This region of the embedded moonlet's gravitational 
influence is of the order of a few Hill radii 

h = a ( — ^ ) , (1) 

\3{m m +m s )J 

where is the semimajor axis of the moonlet, m m its mass and m s the mass 
of Saturn. Compared to the moonlet's semimajor axis the Hill radius is usually 
very small. For large propeller moonlets, like Blcriot or Ear hart, the ratio h* — 



h/do is approximately 10 (Tiscareno et al. 2010). This low ratio naturally 



allows the splitting of the rings into the two regions, where the scattering region 
is shrunk to a line. 

For the rest of the ring, where the moonlet's gravity is negligible, the moonlet 
induced structures are assumed to relax due to inelastic collisional cooling and 



viscous diffusion (Spahn and Sremcevic 2000 Sremcevic, Spahn, and Duschl 



2002|. In previous hydrodynamic models the former process, the relaxation of 
the granular temperature T ~ c 2 /3, has not been considered. This thermal 
relaxation and 3D-scattering are the major themes of this work. 

2.2. Encounter with the moonlet - gravitational scattering 

We describe the encounter of ring particles with the moonlet in a corotating 
frame, rotating about Saturn with the mean motion of the moonlet. The frame's 
origin coincides with the mean orbital location of the moonlet. The x axis points 
radially outward, the y axis points into the azimuthal direction and the z axis is 
normal to the ring plane in a way that the axes form a right handed coordinate 
system. 

With this, the dynamics of ring particles in the corotating frame is given by 
r + 2fJ x r + fl Q x (fl Q x r) = -V$ s - V$ m , (2) 
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where $ s and $ m are the gravitational potentials due to Saturn and the moonlet. 

The ring particles are on orbits with low eccentricity and inclination and 
the mass of the moonlet is very small compared to Saturn's mass m m /m s <C 1. 
In the vicinity of the moonlet the equations of motion of ring particles are well 



approximated by Hill's equations (Hill 1878 Henon and Petitl 1986) 



With the scaled coordinates x = x/h, y = y/h, z = z/h and scaled time 
t' = ft t Hill's equations become 

x = 2y + 3x — 3x/s 3 

y = -2k- 3y/i 3 (3) 
5 = — z — 3z/s 3 , 

where s 2 = x 2 + y 2 + z 2 is the scaled distance to the moonlet and x = dx/dt' . 
These equations are point symmetric about the position of the moonlet (x = 
y = z = 0), and quite comfortable, they do not depend on the moonlet mass 
anymore. All information of the moonlet mass is contained in the scaling length 
h. 

When the ring particles are not in the vicinity of the moonlet, i.e. 1/s 3 — > 
and therefore |V$ S | — > 0, their trajectories are well described by the solutions 
to the homogeneous Hill's equations 

x(t') = a — e cos(t' + f) 

y(t') = C - -at' + Te sm(t' + f) (4) 
z(t') = isin(t' +&) . 

The semimajor axis, eccentricity and inclination are scaled according to 

a-a „ eao . ia 

a = — ' e = lT' l = (5) 

where a spherically symmetric planet ($ s ~ —Gm s /r) is assumed. The phases 
f and oj are the longitude of pericenter and the longitude of the ascending node 
respectively. 
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2.2.1. Test particle integrations 

We integrate the equations of motion numerically for a set of test particles 
using a 5th order embedded Runge-Kutta scheme with adaptive step size control 



(Press et al. 1992). The moonlet and Saturn are assumed to be point masses. 
The integrations start upstream of the moonlet at an azimuthal distance of 
1000 Hill radii to the moonlet and are terminated when the test particle's az- 
imuthal distance exceeds 1000 Hill radii. These limits ensure that the region of 
interaction is well included in the integration. 

We take the initial semimajor axis of the test particles to be uniformly 
distributed in the range of —20 to 20 Hill radii radial distance to the moonlet. 
The initial eccentricities and inclinations of the test particles are choosen to be 
Rayleigh distributed 



f(e, i) = ^ exp 



2c 



^exp 



2c 



r-2 



(6) 



with uniformly distributed initial phases ip and £. This is a fair assumption for 



many kinds of disks (Petit and Henon 1987 Ida and Makino, 1992 Lissauer 



1993 Ohtsuki and Emori 2000). 



This distribution of orbital elements corresponds to a triaxial Gaussian ve- 
locity distribution 



f{Vx,Vy,V z ) 



TIC 



1 

r^exp 



■40.1 



1 



=5? exp 



with a diagonal velocity dispersion tensor T, where T x 



(7) 



fe/2) 2 



and T zz = 



The scaled quantities c Xl c z are related to the unsealed ones by 



c x = c x /(Ct h), c z = c z /(tt h). (8) 
We choose the ratio c z /c x = 0.65, which is consistent with the above Gaus- 



sian velocity distribution (Goldreich and Tremaine 19781 and which is a rea- 



sonable value if we neglect self gravity (Salo et al. 2001 ). 

This deviation from the Maxwellian velocity distribution of conservative 
gases is an expression of the non-equilibrium state, caused by the steady dissi- 
pation of orbital energy in collisions. This energy loss is counterbalanced by the 
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energy gain due to viscous Kcplerian shear, guaranteeing a quasi-equilibrium 
manifested by a steady state. This energy pumping goes at the expense of the 
collective Keplerian motion and thus the ring particles spiral slowly into Saturn 
in the long term (10 s years). 

During the trajectory integrations, the test particles and the moonlet are 
assumed to be point masses and the minimal distance of the test particles to 
the moonlet is recorded. For simplicity, we ignored particles which would collide 
with a spherical moonlet of finite size in further calculations. To estimate the 
difference, we also examined the case where all particles are used. Although 
the radial position of the maximally induced ring temperature changes a bit, we 
found no significant changes in the resulting gap averaged heights. 

2.2.2. Mass transfer 



We use the approach of Spahn and Wiebicke (1989) to describe the mass 



transfer through the scattering region. Motivated by the chaotic behaviour 



of single particle trajectories near the moonlet (Petit and Henon 19861, the 
encounter of ring particles with the moonlet is modelled by a probabilistic 
Markov chain model. The results of the test particle integrations are used 
to calculate transitional probabilities between initial (x, z) and final (S/,z') po- 
sitions of the test particles at the azimuthal boundary of the scattering region. 
These probabilities define a scattering operator A, where A(x' , z' \ x, z) dx dz 
is the probability that matter is scattered from + dx) x (z,z + dz) to 
(x' , x' + dx') x (z' , z' + dz'). The primes denote values after the scattering. 

We assume that the scattering region can be approximated by the x-z plane 
at y = 0, which connects Saturn and the moonlet and is analogous to the 



scattering line of Spahn and Wiebicke ( 1989 ) . In our model, the complete 



action of the moonlet on the ring particles happens at this scattering plane. 
The scattering operator relates the azimuthal mass flux entering the scattering 
region to the azimuthal mass flux leaving the scattering region 



Jjj(x',y = , z')\ = // dx dz A(x', z' | i, z) | Jy(x, z, y — T ) | . (9) 
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Mass conservation is expressed as a condition for the scattering operator 



dx dz A(x , z | x, z) = 1 



(10) 



also illustrating that accreting particles are ignored here. 

For the numerical calculation of the scattering operator we divide the ra- 
dial and the vertical direction into bins. The bin (J, k) describes the region 
(xj, 5j+i) x (zfc, Zk+i), j being an index in radial direction and k one in vertical 
direction. In the following, primed bin indices refer to the situation after the 
scattering by the moonlet and unprimed ones to the situation before. 

The discretized scattering operator is then calculated by 

N(j',k',j,k) 



A(j',k'\j,k) 



(11) 



N(j, k) 

where N(j,k) is the number of test particles starting in bin {j,k). We use 
an averaging procedure to calculate the number of test particles N(J' ,k' ,j,k) 
starting in bin (j,k) and ending in bin (j',k'). For each radial end bin j', 
describing the interval (Xj>, a^'+i), the time per orbit, the test particle will stay 



in this radial end bin is calculated (Spahn 1987) 
At(j') _ 



P(f) 



arccos 



arccos 



N(f , kl , j, k) is then calculated by 



(12) 



(13) 



where the sum is over all test particles and 5^ n '(j,k) is 1 for test particles 
starting in bin (j, k) and otherwise. 

Outside the scattering region we describe the ring using hydrodynamical 
equations. We assume that the fluid parcels entering and leaving the scattering 
region have a Keplerian azimuthal mean velocity Uy(x) — — 3J/2. The azimuthal 



mass flux at the scattering plane is then ( Spahn and Wiebicke| 1989) 
Jy(x, y = ± ,z) = p(x, y = ± , z) u y (x) , 



(14) 
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Figure 1: (a) Mass density p scaled by £o//i after scattering by a moonlet with a Hill radius of 300m. The 
contour lines are for p\ — 1.15 po(z = 0), P2 = Pq{z — c z ) and P3 = po(z = v3c z ), where po, given by equation 



(16), is the mass density before the scattering by the moonlet. Clearly seen is the less dense gap region between 
1.5 < \x\ < 4. (b) Ring thickness after the scattering by the same moonlet measured by: standard deviation a z of 
the vertical mass density profile, root mean square z rms of the vertical excursions and vertical thermal velocity c z , 
all at y = + . All three are within a few percent of the values before the scattering by the moonlet, illustrating 
that the moonlet's gravity alone does not induce considerable vertical excursions of ring particles. 



where p denotes the mass density. The mass transfer through the scattering 
region is thus calculated by 



p(f, k') = A W> k ' I 3, k) p(j, k) \uy{j)\ 

\ u vU )\ jjk 



(15) 



We assume the ring to be in an equilibrium state before the encounter with 
the moonlet. In the following all quantities of the unperturbed ring have the 
subscript 0. The surface mass density So shall be constant, so that the mass 
density in the thin-disk approximation is given by 

E /h 



po(x,y, z) 



exp 



zc z0 



(16) 



The two plots in Figure [T] are made for a moonlet with 300m Hill radius and 



the orbital parameters of Earhart (Tiscareno et al. 20101. We have choosen the 
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velocity dispersion of the unperturbed ring to be cq = 3.9mm/s, corresponding 
to c z = 0.052. Because of the point symmetry of the Hill equations pjl), we show 
only values for x < in the plots. 



Figure la shows the mass density p scaled by So/ft. after the scattering by 
the moonlet. The mass density in the regions \x\ < 1 and |i| > 5 is nearly 
unchanged, whereas the gap region 1.5 < \x\ < 4 has a considerably lower 
mass density than the ring before the scattering by the moonlet. There are 
two regions of high mass density around = 1 and \x\ = 4.25. The solid line 



= 0.0) in Figure 3b shows the corresponding surface mass density. 

Also shown are three contour lines, corresponding to characteristic values of 



the unperturbed mass density po given by equation (16). The first contour line 
is for pi = 1.15 • po(z = 0), i.e. 15 percent larger than the maximal value of pq. 
The second contour line is for the mass density value pi = po(z = c z q). For low 
optical depths c z is a good estimate of 5 rms , the root mean square value of the 
vertical excursion of the ring particles. The third contour line shows the mass 
density p$ = po(z = \/3c z o), where z = \[^>c z is half of the effective geometric 
thickness of the ring. The effective geometric thickness H e g = Vl~2z rms w \fvic z 
of the ring corresponds to the width of a uniform vertical mass density profile 
with the same standard deviation as a Gaussian one. 



Figure lb shows different measures of the ring thickness after the scattering 
by the moonlet. To calculate the ring thickness from the mass density, we 
consider the mass in bin (_?', k') divided by the total mass in bins with bin index 

where we assumed uniform bin sizes. As a quantity, q describes the distribution 
of mass in the vertical direction and formally behaves like a probability, i.e. 
< q(j',k') < 1 and q(j' ,k') — 1. The vertical displacement of the ring 
plane from z = 0, calculated from the mass density, is then given by 



Pz,j> 



I>'^"' fc '), (18) 
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with Zk' being the midpoint of the vertical interval (zk 1 , Zk'+x)- The scale height 
of the ring is then proportional to 

*»(j')= (Yl(^-fi z ,r) 2 q(f,k')\ . (19) 

We also calculated the root mean square value of the vertical excursion of 
ring particles directly from the results of test particle integrations. With the 



averaging procedure of equation (12 1 z rms is given by 



with 

,~ x _E„* (n) p (n) (i') 



1/2 



Zrms,j' - ^ „fnW«V\ ' ^ U > 



(21) 



For the unperturbed mass density po, both a Z Q and z rmsQ equal c z q. The val 



ues after the scattering by the moonlet shown in Figure lb are within 6 percent 
of the values of the unperturbed ring. Calculations for moonlets with different 
Hill radii from 50m to 500m confirm that the difference of the unperturbed 
values to the ones after the scattering by the moonlet are small, e.g. a few 
percent. 

Up to now we have only considered the gravitational interaction with the 
moonlet. Thus, the values of a z , determined from the mass density p, and of 
z rms , the root mean square of the test particle's vertical excursions, are the 
result of gravitational excitement only. The increase of a few percent of these 
values can not explain the height of propeller structures as infered from their 
shadow length. A picture which changes drastically if a moonlet on an inclined 
orbit is considered, which is an issue of ongoing future work. 

Because the moonlet does not induce considerable vertical excursions of ring 
particles, we restrict, in the following, our description to the ring plane and 
model the vertical propeller structure by the granular ring temperature T = 
(Z c + ~cl + ~cl)/Z. 
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Figure 2: Moonlet induced thermal velocities as functions of the radial coordinate x at y = + . All thermal 
velocities are scaled by £l h. (a) Thermal velocities after the gravitational scattering by a moonlet with 300m 
Hill radius. The moonlet induces much more thermal motion in the ring plane than in vertical direction, (b) 
Comparison of the radial thermal velocity c x after the scattering by moonlets with different Hill radii. In the 
gap region —4 < x < —1.5 the radial thermal velocity scales well with Q h for the larger moonlets (h = 100m 
and above). 
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2.2.3. Velocity dispersion 

In this subsection we determine values of the moonlet induced thermal veloc- 
ities after the scattering by the moonlet, which are later used as initial conditions 
for the hydrodynamical equations. We assume that after the scattering by the 
moonlet, the particles are on Keplerian orbits which are well described by equa- 
tion Q. We calculate the thermal speed components c x , c y and c z by taking 
the weighted sample standard deviation of v x , v y — ix/2 and v z for the set of 
integrated test particles. For c x for example 

where f denotes the radial bin number. 



Due to the averaging process ( 12 ) our model does not describe moonlet 



wakes, which are moonlet induced coherent motions of ring particles (Showalter 



et al. 1986 Spahn, Scholl and Hcrtzsch 1994). Nevertheless our model describes 



the gap region quite well, because test particles ending in this region are often 



on trajectories which are very sensitive to initial conditions (Hcnon and Petit 



1986), and their phases arc mixed up. Furthermore, wake structures start to 



dissolve when nearby streamlines cross. For the gap region this happens very 
fast, i.e. for the middle of the gap at |x| = 2.5 after 0.4 orbits. 

For the region beyond the gap however, we have to keep in mind that we 



might overestimate c x and c y with equation ( 22 ) . The averaging process destroys 
the coherent phase relations and systematic motion looks partly like thermal 
motion. In this work we focus on the moonlet induced gap, which is important 
on its own, because of the large thermal excitation in the gap region (cf. Figure 

n. 



The values of the mean thermal velocities, calculated with equation (22), 
are the result of the gravitational interaction with the moonlet only. Figure 



lb shows c z and z rms , which do not deviate much at y = . In Figure 2a 



we compare the mean thermal velocities after the scattering by the moonlet of 
Figure [I] The maximum excitement in each case lies near the inner edge of the 



13 



moonlet induced gap at a radial position of about \x\ = 7/4. 

The moonlet induces much more thermal motion in the ring plane than 
in vertical direction. This can be seen in the ratio c z /c x after the scattering, 
which is several times smaller than the equilibrium value of 0.65. The ratio of 
the horizontal components c v /c x , however, is close to the equilibrium value of 
0.5. Therefore, although the graviational interaction with the moonlet induces 
thermal motion in vertical direction, it is negligible compared to the thermal 
motion induced in the ring plane for moonlets on non-inclined orbits. 



Figure 2b shows the x-component of the velocity dispersion, scaled by flah, 
at y = + after the scattering by moonlets with Hill radii of 50m to 500m. For 
the larger moonlets, with a Hill radius above 100m, the x-component of the 
velocity dispersion in the gap region, especially the maximal value of c x , scales 
well with the Hill radius of the moonlets. 

2.3. Hydrodynamic Flow 

We describe the ring outside of the scattering region with hydrodynamical 
equations. Because of the vertically ineffective gravitational scattering of ring 
particles by the moonlet, we describe the thickness of the ring through the 
ring temperature, and use vertically averaged balance equations. The mass and 
momentum balance are given by 
BY 

— +V-(Eu) = (23) 

^ + (u.V)u = f-ivoP. (24) 

Here E, u, f and P are the surface mass density, the mean velocity, external vol- 
ume forces, and the pressure tensor. The pressure tensor is given in Newtonian 
form 

P = pi 2£idD - SC(V ■ u)I , (25) 

where I is the unit tensor. The vertically integrated pressure, kinematic bulk 
and shear viscosity are denoted by p, £ and v. The traceless shear tensor D is 
given by 

D = - (Vou + uo V) - -(V • u)I , (26) 



14 



where (V o u)jf. — Vjttfc and (uo V)jfc = VkUj. 

The energy balance of the ring particle's random motion reads 



-P : e - V - Q-T 



(27) 



where T = (c 2 , + c 2 + c|)/3 is the granular ring temperature and Q = — K£>VT 
the heat flow in the ring with heat conductivity up. The friction term P : e 
describes the viscous heating of the ring, and T accounts for the energy loss due 
to inelastic collisions. 



2.3.1. Mass and momentum balance 

After the gravitational scattering by the moonlet has opened a gap, the dif- 
fusion of the ring particles, described by the nonlinear viscous diffusion equation 



3d_ 
r dr 



d 



= 



(28) 



for the surface mass density E, will smooth out the induced structures (Spahn 



and Sremcevic 2000 Sremcevic, Spahn, and Duschl 2002) 



In order to solve equation ( 28 1 we apply the following simplifications, already 



used by Sremcevic, Spahn, and Duschl (2002): Let <j\ = S — So, where E is 



the equilibrium value of the surface mass density of the unperturbed ring. In 
the special case of constant kinematic viscosity vq and without curvature terms, 
which is consistent with the Hill approximation of the gravitational scattering, 



equation ( 28 1 reduces to a linear partial differential equation 

d 2 a l 



dip 



where the constant K is defined by 



K 



dx 2 



(29) 



(30) 



(x = y 



Equation ( 29 ) is point symmetric with respect to the position of the moonlet 
z = 0). For the region defined by x < and tp > 0, with boundary 
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Figure 3: (a) Contour plot of the surface mass density downstream of the moonlet. (b) Surface mass density as a 
function of x for different longitudes, illustration the gap-closing. The solid line shows the surface mass density 
after the scattering by the moonlet. 

conditions 



oi (5 < 0, (p = + ) = <Ti4 

<Tl{x — > — OO, ip) = 
01 — > OO) = , 



(31) 
(32) 
(33) 



Sremcevic, Spahn, and Duschl (2002) derived several Greens functions for dif- 
ferent constraints at x = 0, <p > 0. The Greens function, which matched their 
numerical solution best, is an equally weighted superposition of Greens functions 
for the two cases 



(T\(x = 0, cp) = 




We use this Greens function, given by 

G(x, cj>) = %~Zo) (30) 2/3 exp (^^) Bi ((3<£) 



(34) 
(35) 

(36) 
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with the scaled azimuth <p = <p/K, and the Airy function Bi(z), to calculate the 
surface mass density via 

(Ti(5,0)= / dx o cr 1)ini ($ ) G(x, 4>; x ) . (37) 

J —OO 

As initial surface mass density c i } mi we use the vertically integrated mass density 



calculated with the scattering operator (15 1, shown in Figure 3b as solid line 



and numerically integrate equation (37). The two plots in Figure [3] are again 
made for a moonlet with 300m Hill radius and the orbital parameters of Earhart. 
For the value of the kinematic viscosity in Saturn's A ring we used vq = 0.01 



m 2 /s (Tiscareno et al. 2007) 



Figure [3a] shows the surface mass density after the scattering by the moonlet 
as function of x and (p. White contour lines represent surface mass density values 
below So, showing the evolution of the gap. Black contour lines enclose regions 
of enhanced surface mass density compared to So- To illustrate the gap-closing, 
we plotted the surface mass density as a function of the radial coordinate x for 
different longitudes <f>, shown in Figure |3b| 

The azimuthal extent of the surface mass density scales with a$K. Because 
K oc h 3 , the importance of the mass diffusion process, in the first few orbits 
after the scattering by the moonlet, depends strongly on the Hill radius of the 
moonlet. For large moonlets, e.g. h = 500m, the surface mass density stays 
close to the initial value at </> = 0. The smaller the moonlet the larger the 
influence of the mass diffusion in the first few orbits. 

2.3.2. Energy balance 

In this subsection we will consider the relaxation of the moonlet induced 
thermal heating by dissipative particle interaction. The vertically ineffective 
gravitational interaction of ring particles with the moonlet can not explain the 
vertical height of propeller structures. The moonlet induces thermal motion 
mainly in lateral direction. 

But due to collisions between ring particles, the thermal motion induced in 
lateral direction will be converted to vertical thermal motion till the ratio c z /c x 
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reaches its local equilibrium value. This asymptotic value of c z /c x is established 



very fast, after a few collisions per particle ( Hameen- Anttila and Lukkari 1980 



Hameen-Anttila and Salo 1993). We therefore use the granular temperature T 



to model the evolution of the vertical propeller structure, infering the thickness 
of the ring from T and the equilibrium value of c z /c x . 

The thermal energy balance equation in the corotating frame reads 

3 f Fit bt bt ~\ 

- E {^ + (n-no)^ + « r ^} = -P:e-V.Q-r . (38) 

The temperature relaxation of the ring is much slower than the relaxation 
of the ratio c z /c z to its local equilibrium value. On the other hand, the tem- 
perature relaxation is fast compared to the mass diffusion timescale. Therefore 
we consider, for simplicity, only viscous heating and energy loss due to inelastic 
collisions. Furthermore, we assume stationarity to obtain 

3 BT 9 

_(O_n )— = -vCll-k 3 Q T(l-e 2 )T , (39) 

where we used P : e = — ^vVL^Yj/A for the friction term, regarding just the 
dominant Kepler shear. For the cooling term we used T = k^iloT^ — e 2 )ST 



with a constant coefficient of restitution e = 0.5 and with = 1.5 (Stewart et 



al. 1984), corresponding to an energy loss rate 

4oii = -w c /6(l-£ 2 )£c 2 , (40) 

with collison frequency w c = 3f2r, and temperature c 2 = 3T. 

We assume the viscosity to be constant v = v$ on timescales which are 
relevant for the exponential decay of the temperature. With a constant viscosity 
and a constant coefficient of restitution, viscous heating and collisional cooling 
balance at a local equilibrium temperature 

T M - 9 "° n ° (41) 
cq[ ' ~ 4A: 3 (1 - e 2 )r ' [ > 

which depends on the local optical depth, or assuming equally sized ring particles 
with particle radius R p and mass m p , via r = 7ri? 2 E/r7i p , on the surface mass 
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density S. The temperature of the unperturbed ring before the encounter with 
the moonlet is Tq = T oq (To). 

The temperature independent viscosity seems to be a strong restriction and 
is considerd to be a first step towards a more complete description of the prob- 
lem. The applicability of our model presented in Section [4] will justify the 
simplifications applied in this work. 

With these assumptions equation ( |39| ) becomes a linear ordinary differential 
equation and with (ft(r) — £Iq) d/dcp — — 3/2ttox d/dy and T\ = T — Tq: 

x ~ = v n ^ + h 3 r(l - £ 2 )7\ . (42) 
ay 2jo 9 

The solution of this equation can be written in the form 

Ti(5, y) = T 1Mi (x) + F f{y")e^^ dy", (43) 



o 



with the functions 



r - -v 4fc 3 (l - s 2 )tq f* S(gjQ ; 
\x\ S 

The plots in Figure [4] show the downstream temperature relaxation. Again 
they are made for a moonlet with 300m Hill radius and the orbital parameters 
of Earhart, using i/q = 0.01 m 2 /s as value of the kinematic viscosity in Saturn's 
A ring. Figure [4a] shows a contour plot of the ring temperature T scaled by the 
unperturbed ring temperature Tq, illustrating the azimuthal temperature decay 
for different radial regions. Ring particles with \x\ < 0.75 are on horseshoe 
orbits and their minimal distance to the moonlet is large, so that there is no 
considerable thermal excitation. The gap region, on the other hand, is highly 



excited for quite an azimuthally extended range. Figure 4b shows the azimuthal 
temperature evolution, at fixed radial position x — —2.0, downstream from the 
moonlet. Also, the local equilibrium temperature T eq is shown as function of the 
surface mass density. For the first few orbits, the ring temperature decreases 
exponentially to the local equilibrium temperature. Afterwards, the viscous 
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xth orbits 

(a) (b) 

Figure 4: Downstream temperature relaxation for a moonlet with a 300m Hill radius, (a) Contour plot of the 
temperature decrease downstream of the moonlet. (b) Temperature relaxation at radial position x — —2. The 
solid line shows the temperature relaxation of the ring, the dashed line corresponds to the equilibrium temperature 
T oq , for which viscous heating and collisional cooling are in balance. 
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1000 
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Figure 5: (a) Decay function 7 at fixed radial position x = —2 for two moonlets, one with 100m Hill radius and 



one with 500m Hill radius. The solid lines show 7 calculated according to equation (44). In contrast the dashed 



lines are calculated with an azimuthally constant surface mass density £ equal to Ei n i- (b) It is shown how many 
orbits it takes the temperature to decay to 5 percent of T eq for moonlets with different Hill radii at x = —2. 
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heating and the collisional cooling are balanced and the ring temperature evolves 
as a function of the surface mass density of the ring. In this case the temperature 
gets within a margin of 5 percent of T eq in about 7 orbits. 

The influence of the mass diffusion on the temperature decay depends on 



the Hill radius of the moonlet. Figure 5a shows decay functions 7, calculated at 
radial position x = —2, for two moonlets with Hill radii of 100m and 500m. The 



solid lines show 7 calculated according to equation (44 1. In contrast the dashed 
lines are calculated with an azimuthally constant surface mass density X equal 
to Sini- For the larger moonlet the solid and the dashed line stay fairly close 
to each other, because the change of the surface mass density is rather small 
on the timescale of the exponential decay. For the smaller moonlet both lines 
differ clearly indicating that the influence of mass diffusion is by far larger. 
A fair approximation to the solution ( 43 1 is given by 

T approx (5;, y) = (T ini - T oq (r)) e~<^ + T oq (r) . (45) 



For > 1 the maximal difference of equation (43) to the above approximate 
solution is about 5 percent, tested for moonlets with Hill radii of 50m to 500m. 
The exponential decay of the temperature stops after a few orbits, when 



the local equilibrium temperature T cq is reached. Figure 5b sketches how many 
orbits it takes to get within 5 percent of T cq for moonlets with different Hill 
radii at fixed radial position x = — 2. Using the approximate solution T applox , 
the azimuthal position y^ where the the ratio T/T cq equals 1 + £ is 

1 ( Ti„i - T cq (r) 



j -s h lwJ ' (46) 

where we defined formally j(x,y) = (7)2/4- Both (7) and T cq depend on the 
optical depth r and therefore on the azimuthal position y. For large moon- 
lets, close in size to ring-moons able to open a circumferential gap, r is nearly 
constant on the timescale of the exponential temperature decay. Because of 
21ni ^> T C q and Tj„i oc h 2 for large moonlets, y^ will scale with y^ oc 21n(/i). 



In Figure 5b a fit of the function y^ = A ]n(Bh ) to the values of the large 



moonlets with Hill radii of 200m to 500m is presented in order to underline this 
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(a) 



(b) 



Figure 6: Azimuthal relaxation of the propeller height for a moonlet with 300m Hill radius, (a) Propeller heights 
as function of radial coordinate x for different longitudes, (b) Downstream relaxation of the propeller height. 
The solid line shows the relaxation at x = —2. The dashed line shows the propeller height radially averaged over 
the gap region. 



scaling. 



3. The vertical height of propellers 

In the thin-disk approximation and with a z independent vertical velocity 
dispersion c z , the vertical profile of the mass density can be described well by 
a Gaussian with standard deviation c z /Qq ) |Stewart et al. 1984 Simon and 



Jenkins 1994 Schmidt et al. 1999). We use the effective geometric thickness 



H e g as a measure of the ring thickness and determine the height of the propeller 
structure as half the effective geometric thickness of the ring. 

To determine c z as a function of the temperature T, we use the equilibrium 
values of the ratios c z /c x and c y /c x . A perturbed ring will reach the equilibrium 



values of these quantities after a few collisions per particle (Hameen-Anttila and 



Lukkari 1980 Hameen-Anttila and Salo 1993). In Saturn's A ring with r = 0.5 



and with a collision frequency of u c = 3f2r, there are about 10 collisions per 
particle per orbit. We assume that after 5 collisions, corresponding to half an 
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h/m 


50 


100 


200 


300 


400 


500 


H eap /h 


0.90 
0.80 


0.90 
0.72 


0.94 
0.72 


0.95 
0.73 


0.95 
0.73 


0.96 
0.74 



Table 1: Heights for moonlets with different Hill radii from 50m to 500m. i? max 
is the height at x — —2.0 and y = 10. H gap is the height averaged over the gap 
region at azimuthal position y = 12. 



orbit, the equilibrium values of these ratios are established. With the relaxed 
ratios {c z /c x ) cq = 0.65 and (c y /c x ) cq — 0.5, the z component of the velocity 
dispersion in terms of T is given by 



Cz[T) (cj cq V 1 + (c,/c.) c 2 q + {c z /c x )l q ■ (4?) 

The two plots in Figure [6] are made using the same parameters as were taken 
for Figure |4j i.e. a moonlet with 300m Hill radius, the orbital parameters of 



Earhart and a kinematic viscosity of uq = 0.01 m 2 /s. In Figure 6a the height 
of the propeller structure as a function of the radial coordinate x is plotted 
for different values of y. The largest heights are reached in the region between 
x = —4 and x = —1, the maximum being at about x — —2 for y = 10 and 
slowly moving to x = —2.5 at y = 100. 

The azimuthal relaxation of the height of the propeller structure, for x = —2, 



is shown in Figure 6b We assume that after half an orbit or 5 collisions per 
particle the equilibrium value of c z /c x is established and that afterwards the 
ring temperature describes the propeller height well. At radial position x = —2, 
half an orbit corresponds approximately to y = 10. The height at y = 10 is 0.95 
Hill radii or about 285m. 

Furthermore, Figure [6b| sketches the propeller height radially averaged over 
the gap region. We choose the approximate middle of the gap at x = —2.5 as 
a reference, and in that case half an orbit orbit corresponds to y = 12. At that 
position the gap averaged propeller height is 0.72 Hill radii or about 215m. 
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Tabic [T] shows propeller heights for moonlets with different Hill radii. H max 
is the height at radial position x = — 2 and azimuthal position y = 10. H gap 
is the height averaged over the gap region at azimuthal position y = 12. The 
maximal heights are close to one Hill radius for the tested moonlet sizes. The 
values of the gap averaged height for the moonlets with Hill radius 100m and 
above are approximatly 0.7 Hill radii. The averaged propeller height for the 
moonlet with a 50m Hill radius is larger (0.8h), which is due to the large value 
of the unperturbed height of the ring scaled by the Hill radius of the moonlet. 



4. Discussion 

We study the vertical extent of propeller structures in Saturn's rings, focus- 
ing on the propeller gap region. The effective geometric thickness H c g is used to 
describe the propeller height as a function of the ring temperature, calculated 
using the equilibrium ratios of the thermal velocities. A vertically constant ring 
temperature serves as a fair assumption, because for the low optical depths 
in the propeller gap region the vertical dependence of the ring temperature is 



rather weak (Schmidt et al. 19991 



In our model the azimuthal temperature decrease is caused dominantly by 
the disturbed balance of viscous heating and collisional cooling. We assume a 
constant coefficient of restitution for the cooling term and constant viscosity for 
the viscous heating term. This is a simplification, but allows a semi-analytical 
solution of the azimuthal temperature relaxation (some integrals have to be 
numerically evaluated), and most important, it is consistent with images of 
propeller shadows taken by the Cassini spacecraft. 

The restriction of a constant vq could be dropped by introducing e.g. a 



power law dependence (Spahn et all 2000) 



U) UJ UJ ' ,48) 

For small moonlets (h < 50m) one has a rather small ratio T^/Tq < 5. The 



power law dependence (48 ) with a — 1/2 (Salo et al. 2001 ) then gives a viscosity 
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about twice as large as the unperturbed one. On the other hand, the surface 
mass density in the gap region after the scattering by the moonlet is about half 
of the unperturbed value, so that a constant viscosity vq can be justified. In 
case of large moonlets the ratio Tj n i/To is by far higher (190 for a moonlet with 
300m Hill radius at \x\ = 2). The viscosity will rise considerably with increasing 
temperature, tending to an increased heating, and thus, to an increased cooling 
time scale. However, because of an increasing coefficient of restitution and due to 
increased collision rates, an enhanced cooling works against that effect. Insofar, 
as a first step, we have used a constant viscosity and a constant coefficient of 
restitution, interpreting both as effective values. 

Furthermore, we have neglected the temperature decrease due to heat con- 
duction, based on the grounds that the relaxation of the ring temperature by 
granular cooling is the faster process, compared to the large diffusive time scales 
blurring the propeller gaps. However, heat conduction as well as viscous dif- 
fusion could become important for small moonlets. We have estimated the 
temperature decrease due to heat conduction for a simplified model, described 
in Appendix A For small moonlets (h < 50m) the temperature decrease due 



to heat conduction is comparable to the temperature decrease of our granular 



cooling driven solution (43 1. On the other hand, for large moonlets the tem- 
perature decrease due to heat conduction was approximatly 20 times smaller 
than our cooling driven relaxation after 10 orbits, justifying a neglect of heat 
conduction for large moonlets, and thus, our approach. 

In the summer of 2009, near Saturn's equinox, the Cassini spacecraft took 
images, which show prominent shadows cast by propeller moonlets. The height 



of the propeller features were calculated from the observed shadow length ( Tis- 



careno et al. 2010 ). For the propeller moonlet Blcriot, a height of 430±30m was 



determined, for Santos-Dumont 120m, and for Earhart 260m. Figure [7] shows 
the propeller moonlet Earhart orbiting near the Encke gap, where Figure |7b| 
portraits Earhart in August 2009, a few days after Saturn's equinox, to cast a 



350km long shadow. In contrast, Figure 7a shows Earhart in April 2008 long 
before equinox. 
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(a) 



(b) 



Figure 7: The propeller moonlet Earhart near the Enke gap. (a) This image was taken by the Cassini spacecraft 
on April 11, 2008 with Cassini's narrow angle camera. The resolution is about 2 km per pixel. The propeller 
structure is about 5 km in radial dimension and about 60 km in azimuthal dimension, (b) Image of Earhart near 
Saturn's equinox, taken by Cassini's narrow angle camera on August 13, 2009. The resolution is 7 km/pixel. The 



shadow cast by the propeller is 350km long and the height of the propeller is estimated to be 260m (Tiscareno 



et al. 2010). Credit: N AS A/ JPL/Space Science Institute 
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Using our values of the gap averaged propeller height from Table [Tj the 
height of the propeller structure is approximately 0.7 times the Hill radius of 
the propeller moonlet. This gives 370m for the Hill radius of Earhart and 615m 
for the Hill radius of Bleriot. Assuming a sperical shape, the moonlet radius 
of Earhart is then 240m for a mass density of 600kg/m 3 , and 280m for a mass 
density of 400kg/m 3 . 

We determined the azimuthal extent of Earhart from Figure [7b] to be 120 ± 
20km, the Cassini ISS team reported about 130km on their website. The pro- 
peller structure thus azimuthally extends about 60 ± 10km downstream from 



the moonlet. In Section 2.3.2 we used a ring temperature of 1.05 • T eq to deter- 
mine the end of the exponential temperature decay. Using this criterium for a 
moonlet with 370m Hill radius and for the radial ring particle position \x\ = 2.5 
(approximately the middle of the gap region) , the exponential temperature de- 
cay stops after about 7 orbits. These 7 orbits at |i| = 2.5 corresponds to an 
azimuthal propeller length of 61km downstream of the moonlet which agrees 



astonishingly well with the 60 ± 10km determined from Figure 7b To calcu 



late the mean decay constant (7), we used equation (46) with the determined 
azimuthal propeller extent, to find (7) = 0.05. At \x\ — 2.5 this corresponds to 
fc3(l — e 2 ) = 1.1, which nicely compares to ^3(1 — e 2 ) = 1.125 used in our model. 

5. Summary and Outlook 

In the present work we study the vertical extent of propeller structures in Sat- 
urn's rings. Our focus lies on the gap region of the propeller and on non-inclined 
propeller moonlets. In order to describe the vertical structure of propellers we 



extend the model of Spahn and Sremcevic ( 2000 1 to include the vertical direc 



tion. We use the scattering operator concept (Spahn and Wiebicke 19891 to 
model the gravitational interaction of the embedded moonlet with the ring par- 
ticles, taking place in a relatively small region around the moonlet. Outside this 
region we describe the Keplerian ring flow using hydrodynamic equations. The 
azimuthal relaxation of the propeller height is modelled by the disturbed bal- 
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ance of viscous heating and granular cooling using the energy balance equation 
of the ring particle's thermal motion. 
Our conclusions and results are: 



1. Considering only the gravitational interaction of ring particles with a non- 
inclined moonlet, the moonlet does not induce considerable vertical excur- 
sions of ring particles. However, the moonlet causes considerable thermal 
motion in the ring plane. The maximal value of this thermal motion, when 
measured by the radial thermal velocity, scales well with the Hill radius 
according to c x ~ hVL. 

2. Collisions between ring particles are an effective mechanism to locally 
convert the lateral induced thermal motion into vertical excursions of ring 
particles. After the scattering by the moonlet, the ratio of vertical to 
radial thermal velocity is far from its equilibrium value, but the local 
relaxation to the equilibrium value takes only a few collisions per parti- 



cle (Hameen-Anttila and Lukkari 1980 Hameen-Anttila and Salo 1993). 



Consequently, after a few collisions per particle, the ring temperature de- 
scribes the vertical extent of the propeller structure astonishingly well. 

3. We use half of the effective geometric thickness of the ring as a measure 
for the propeller height in the gap region and find maximal heights on the 
order of the Hill radius of the moonlet. The height, averaged over the gap 
region, is obtained to be about 0.7 Hill radii. 

4. For the first few orbits, the azimuthal evolution of the ring temperature is 
characterised by an exponential decay to a local equilibrium temperature. 
Later, the viscous heating and the collisional cooling are balanced and the 
ring temperature evolves as a function of the local optical depth of the ring. 
The fast decay is consistent with images of Earhart, taken by the Cassini 
spacecraft near Saturn's equinox. We found that for the exponential decay 
phase the neglect of heat conduction is a useful and valid assumption for 
large propeller moonlets like Bleriot and Earhart. For small moonlets, e.g. 
with Hill radius of 50m, heat conduction and dissipative effects become 
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important right after the scattering. 

Future work should incorporate a kinetic description of the counteracting 
viscous heating and collisional cooling, to study the consequences of a velocity 
dependent coefficient of restitution and a temperature dependent viscosity on 
the azimuthal relaxation of the ring temperature. 

Another subject of ongoing work is the vertical extent of the wake region of 
propellers. By dropping the averaging process in the calculation of the scattering 
operator, radial and vertical phase information of the ring particle's motion is 
preserved. This allows us to describe the coherent motion of ring particles in the 
wake region. An extension to inclined moonlets would then allow comparison 
to the wake region of the ring-moon Daphnis. 

Further work should be done on the inclusion of heat conduction, viscous 
diffusion and fluctuations into our model, to better describe propellers of small 
moonlets. Finally, we plan to compare our results to direct N-body simulations 



of propellers (Seifi et al. 2005 Sremcevic et al. 2007). 
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Appendix A. Influence of heat conduction on the temperature decay 



We neglected the heat conduction term in equation ( 39 ) on the grounds that 



the relaxation of the ring temperature is, compared to the mass diffusion, a fast 
process. To check for which parameters this assumption is valid, we use our 



solution of (|39|) to estimate how large a heat conduction term in equation (|42 1 
i have 



would have been. With a heat conduction term, equation (42) reads 

4k (d 2 T x d 2 T x 



<j\ 4 
^— - vqVLq — + -A: 3 t(1 
dy E 9 



9fi h 2 \ dx 



dy 2 



(A.l) 
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with the kinematic heat conductivity re, which we assumed to be constant. Fur- 



ther we assume k/v 20 (Salo et al. 20011. The azimuthal part of the heat 
conduction term is several times smaller than the cooling term which describes 
the inelastic collision between ring particles. For a moonlet with 50m Hill ra- 
dius, it is at least 20 times smaller, whereas it is at least two orders of magnitude 
smaller for a moonlet with 200m Hill radius. The radial part of the heat con- 
duction term, on the other hand, is of the order of the cooling term for |x| < 1.5. 
For small moonlets, e.g. h = 50m, this applies also for |x| < 4. 

On the timescale of the mass diffusion process, heat conduction can not be 
neglected, but for small moonlets heat conduction seems to be also important 
in the first few orbits. Therefore, we estimated the temperature decrease in the 
middle of the gap at x = —2.5. For that, we assumed a Gaussian radial tem- 
perature profile at y = + , centered at xq = —2.5 and with standard deviation 
a calculated according to full width at half maximum of 3: 

T(x, 0) = r(3o, 0) exp (- {x "j^H . (A.2) 



Using this as initial condition, wc solved equation (A.ll neglecting all but the 
heat conduction term 

X ° dy 9fl h 2 \ Dx 2 dy 2 J 1 [ ' 

further simplified by setting x = xq. The azimuthal evolution of the temperature 
at radial position Jo ist then given by 

it; \ -V 2 



T( X0 ,y) = V^T M {2^ + ^^y) . (A.4) 

For a moonlet with a Hill radius of 50m, the temperature decrease according to 
the above solution was approximatly equal to the temperature decrease of our 



solution (43 1 in the first few orbits. For a moonlet with 200m Hill radius on the 
other hand, the temperature decrease due to heat conduction after 10 orbits was 
approximatly 20 times smaller than the temperature decrease of our solution 



(43). This indicates, although we neglected the Kepler shear by setting x = xq 



in equation (A.3), the importance of heat conduction for the temperature decay 
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for small moonlets. For large moonlets, 200m Hill radius and above, the neglect 
of heat conduction seems to be a valid assumption. 
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